Covariant cosmological perturbation dynamics in the large-scale limit 
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Using the existence of a covariant conserved quantity on large perturbation scales in a spatially 
flat perfect fluid or scalar field universe, we present a general formula for gauge-invariantly defined 
comoving energy density perturbations which encodes the entire linear perturbation dynamics in a 
closed time integral. On this basis we discuss perturbation modes in different cosmological epochs. 
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I. INTRODUCTION 

Because of its notable significance for the early stages 
of structure formation in the universe cosmological per- 
turbation theory has received considerable attention dur- 
ing the past decades (see, e.g. 

Depending on the choice of basic perturbation vari- 
ables there are different possibilities to formulate the 
dynamics of primordial inhomogeneities. These choices 
may refer to different gauges or to different sets of gauge- 
invariant variables. In this paper we focus on the covari- 
ant approach which uses exactly defined tensorial quan- 
tities ( instead of conventional, generally gauge- 
dependend perturbation variables. A covariant approach 
is conceptionally superior to noncovariant approaches 
since it avoids the explicit introduction of a fictitious 
background universe and thus circumvents the gauge- 
problem, which just originates from the nonuniqueness 
of the conventional splitting of the spacetime into a ho- 
mogeneous and isotropic zeroth order and first-order per- 
turbations about this background. The basic dynamic 
quantity used by Olson ||, Woszczyna and Kulak 
Ellis and Bruni M and Jackson Q is the covariantly de- 
fined spatial gradient of the energy density of a comoving 
(with the four-velocity of the cosmic fluid) observer. The 
corresponding dimensionless, fractional quantity obeys a 
second-order differential equation (see, e.g. 0). 

As was shown recently |9| , the linear cosmological per- 
turbation theory of an almost homogeneous and isotropic 
universe may be simplified by introducing covariant vari- 
ables defined with respect to hypersurfaces of constant 
expansion, constant curvature or constant energy den- 
sity. These new quantities which represent suitable lin- 
ear combinations of Ellis-Bruni-Jackson type variables 
turned out to be particularly useful to characterize a con- 
served quantity on large perturbation scales which was 
applied subsequently to study the perturbation dynam- 
ics in an inflationary universe. pQ|. Generally, the im- 



portance of conserved quantities is well recognized in the 
literature ( pd|-p^]). In this paper we clarify how the 
existence of a covariant conserved quantity may be used 
to solve the perturbation equations in the corresponding 
limit. Our main purpose is to derive a general formula 
which reduces the entire large-scale linear perturbation 
dynamics for arbitrary equations of state to a closed time 
integral over homogeneous "background" quantities. We 
then demonstrate how the cosmological mode structure 
for all cases of interest (exponential inflation, power-law 
inflation, radiation and matter dominated periods etc.) 
follow from the general formula in an elementary way. 

This paper is organized as follows. Section II briefly 
recalls the basic relations of a covariant perturbation the- 
ory with special emphasis on the characterization of a 
conserved quantity in the large-scale limit. Section III is 
devoted to the general solution of the large-scale pertur- 
bation dynamics while section IV presents the detailed 
mode structure for a variety of cosmologically relevant 
equations of state. A short summary is given in section 
V. 



II. BASIC RELATIONS OF A COVARIANT 
PERTURBATION THEORY 

We consider the cosmic medium characterized by an 
energy-momentum tensor with perfect fluid structure, 



(m,n... = 0,1,2,3) , (1) 



where p is the energy density, p is the pressure, u m is the 
four- velocity {u m u m = -1), and h mn = g mn + u m u n the 
spatial projector on surfaces orthogonal to the 4- velocity. 
We recall that the structure (Q) is also valid for a mini- 
mally coupled scalar field with the identifications 



p = \^-V{4>) 



(2) 
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and 



where <fi = 4>, a u a = \/ —g ab 4>,a4>.b and V (<f>) is the scalar 
field potential. 

Local energy momentum conservation T zk k = 0, implies 



The three-curvature in this case is known to be 7^ <0> = 
6k/a 2 , where k = 0, ±1 and a now coincides with the 
(3) scale factor of the Robertson- Walker metric. 

In order to describe the dynamics of the inhomo- 
geneities in terms of D a , P a and t a we differentiate the 
first Eq. (Q), project orthogonal to u a and multiply by 
a. The left-hand side of the resulting equation may be 
written as (cf. (§)) 



p = -e (p + P ) , ( P + P )u m = - Pik h r ' 



(4) = Kn - au m p + (w c m + cr c m ) ah™p. 



Here, = u\ is the fluid expansion and u m = u"^u n the 
fluid acceleration. 

Additionally, we make use of the Raychaudhuri equa- 
tion 

9 + ie 2 + 2 (a 2 - lo 2 ) - u« - A + I (p + 3p) = , 

(5) 

where A is the cosmological constant and k is Einstein's 
gravitational constant. The magnitudes of shear and vor- 
ticity are defined by 



with 



^■2 — ^ ,-r ab , ,2 — , , , ,af> 



1 



(6) 



Cab = h c a hb u (c;d) - -®h ab , w ab = h c a h b U[ c . d ^ . (7) 
The 3-curvature scalar of the projected metric, 

11 = 2 (-^O 2 + a 2 - uj 2 + up + A^j , (8) 

reduces to the 3-curvature of the surfaces orhogonal to u a 
for lu — 0. The homogeneous and isotropic Friedmann- 
Lemaitre-Robertson- Walker (FLRW) universes corre- 
spond to the special case a = to = ii a = 0. 

Suitable covariant variables to characterize spatial in- 
homogeneities are [(7j-[io[ 



D„ 



Pa= a -^, t a ^a K e, c , (9) 
P+P P+P 



where a is a length scale generally defined by O = 3d/ a. 
The quantities D a and P a represent fractional, comoving 
(with the fluid four-velocity) energy density and pres- 
sure perturbations, respectively. Inhomogeneities in the 
expansion are described by the quantity t a . 

In a homogeneous and isotropic FLRW universe (su- 
perscript 0) the Friedmann equation 

(0) 1 /„(°>\ 2 1 _,(0) 

k P = -(e j + -n -a, (io) 

holds and the Raychaudhuri equation (||) reduces to 

A «» 3 / CO) (0)\ 1(0) 

+ 2 V P J = 2 ■ ^ ' 



(12) 

where we have used the well-known decomposition 
- _■ @ h 

"Ui-.n — U%Ufi -\~ (Jin ~T LOi n -\- ^ tlin 

of the covariant derivative of the four- velocity. 

Restricting ourselves to first-order deviations from ho- 
mogeneity and isotropy we find [pi 



D 



P 



D M + t M = 0, {p,v... = 1,2,3) . (13) 



Analogously, one obtains from Eq. (||), 

2„. K , , „ (\„ V 2 
3 



^ = --6^--^+^)^- I )P„ . (14) 



Combining the last two equations and defining the sound 
velocity c s in a standard way as c 2 = p/p, the first-order 
inhomogeneities are governed by the second-order equa- 
tion 



D„ 



(c 2 )-e 



+ U(P- 3 P) + 2A C '+ o (P+P) 



^p 



(15) 



ID 
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For a fluid (superscript f) one has P a 
Eq. ( |l5| ) corresponds to Jackson's || equation (57). 
a scalar field (superscript s), because of h c a 4>, c = 0. 

potential term neither contributes to D a ' nor to P l a 



and 
For 
the 



and, consequently, P^ = D K °' is valid |10| which, if used 



in Eq. (pL 5|) , results in a closed equation for D n as well. 

The well-known second-order differential equation ( |l5| ) 
together with a corresponding expression for P u provides 
us with a complete description of linear scalar perturba- 
tions in perfect fluid and scalar field universes. Our aim 
here is to find the general solution of equation (|l5|) in 
the large-scale limit for spatially flat (fc = 0) cosmologi- 
cal models. To this purpose we rewrite the perturbation 
dynamics in terms of a different basic variable. 



III. SOLVING THE LARGE-SCALE 
PERTURBATION DYNAMICS 

As follows from the definitions (^|), the quantities D a , 
P a and t a are defined with respect to comoving hypersur- 
faces. It was shown in Q that the perturbation dynami- 
cal description simplifies if written in terms of covariant 
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variables defined with respect to hypersurfaces of con- 
stant curvature, constant expansion, or constant energy 
density. These variables correspond to certain combina- 
tions of Ellis-Bruni- Jackson type variables. For example, 
the exactly defined covariant quantity 



£>i ce) = D„ 



P + PO 



(16) 



represents in first order the fractional, spatial gradient of 
the energy density on hypersurfaces of constant expan- 
sion (superscript ce) P,fTo|]- Rewriting the linear pertur- 
bation dynamics in terms of D a ce ^ and restricting our- 
selves to k = 0, one obtains 



(17) 



Except for the spatial pressure gradient on the right-hand 
side of Eq. ( p"7j ) all terms of the dynamical perturba- 
tion equation may be included into a first time deriva- 
tive. With Df, = D( m ) V M Q( m ), P M = P( m )V A( Q( ro ) 

and D^ e ' = d'^V ( m ) (and corresponding relations 
for the other perturbation quantities), where the Q( m ) 
satisfy the Helmholtz equation V 2 Q( m ) = ~ m2 Q(m)i 
the quantity m is related to the physical wavelength by 
A = 2na/m ( |^,[l6|-|l8|] ) . It follows that the spatial gra- 
dient terms on the right-hand side of Eq.(|l7|) may be 
neglected on large perturbation scales (m <§; 1) and the 

quantity a 2 0.D^j is a conserved quantity both for a per- 
fect fluid and a scalar field. 

Denoting the conserved quantity by — -E( m ), i.e. 

a 2 QD[2] = -E {m) = const , (to « 1) , (18) 
and taking into account the relation 



® n(«0 

e " 



between and which follows from the definition 

( |l6| ) of and Eq. (|l3|), the equation to solve is 



e 
e 



(19) 



D 



(m) 



- I - + cfe D {m) = 



v 6 ' -»-]-*•>> a 2 Q ' 
With = 3d /a and the abbreviation 



(m < 1) . (20) 



(21) 



the general solution to Eq. (20) for m <C 1 is 



Z3(m) = Qa 3c °exp{-q} 
t 



dt 



6 2 a 2 + 3c 



exp {<?} + C (r , 



(22) 



where Ct m ) is an integration constant. For many pur- 
poses it is sufficient to consider c 2 w const in which case 
q = and the solution (22) reduces to 



D 



(m) 



6a 



dt 



E, 



(m) 



6 2 a 2 + 3c = 



(m) 



(m < 1) . 
(23) 



Formula (|22| ) comprises the entire large-scale linear per- 
turbation dynamics for arbitrary equations of state. It is 
the main result of this paper. In the following section we 
use the expression ( p2[ ) (or the special case (23)) to derive 
the mode structure for different cosmologically relevant 
equations of state. 



IV. COSMOLOGICAL MODE STRUCTURE IN 
THE LARGE-SCALE LIMIT 

A. Vanishing cosmological constant (A = 0) 

Let us assume equations of state p = (7 — 1) p with 
constant values of 7 in the range < 7 < 2, i.e. c 2 = 
p/p = 7 — 1- The case 7 = 0, equivalent to p = —p 
and c 2 = —1, is characterized by p — const, = 3H = 
const, i.e. a cx exp[Ht\. Use of these dependences in 
formula ( p3|) yields a behaviour oc a~ 2 and oc a~ 3 for 
the cosmological modes. Any perturbation about the de 
Sitter spacetime is exponentially damped JToft - 

For < 7 < 2 we have p oc a~ 37 and, via Eq. (|lC|), 
a oc t 2 ^ 3j \ Since 6 a T 1 one obtains a behaviour 
cx a 3 ^ 2 oc i 2 - 4 /(37) and oc a 37 / 2 " 3 cx i 1 " 2 /? for the 
large-scale modes Q which coincides with corresponding 
results in fli| . 

It follows that in the range < 7 < 2/3 where a oc 
t n with n > 1 (power-law inflation) both perturbation 
modes are decaying. The case 7 = 2/3, corresponding to 
so-called K-matter J|(J, is characterized by a oc t. The 
dominant mode turns out to be constant in this case pl| ] 
while the second mode decays as a~ 2 oc i~ 2 . 

In the interval 2/3 < 7 < 2 there is one decaying and 
one growing mode. For 7=1 (dust) one reproduces the 
well-known behaviour oc a oc t 2 / 3 and oc a 3 / 2 oc t 1 
for the large-scale perturbation modes, while one obtains 



modes oc a 2 oc t and oc a -1 oc i -1 / 2 for 7 = 4/3 (radia- 
tion) . 

For 7 = 2 (stiff matter) there is no decaying mode. The 
first term in (|2^) grows as a 4 oc i 4 / 3 while the second one 
remains constant. 

None of these results is new. Our point here is to 
demonstrate the possibility of a unified representation 
of the cosmological mode structure via formulae ( p2| ) or 
(p3|). We believe this representation to be useful also 
under more general circumstances, where, e.g. c 2 varies 
in time or in multifluid cosmological models to be studied 
elsewhere. 
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B. A-dominated universe 



For A 3> up in Eq. (|1CJ) we have 6 = 3H = const, i.e. 
a oc exp [Hi] . Matter on this background is again char- 
acterized by p — (7 — 1) p with c\ = (7 — 1). ^From 
the general formula (g3j) we find perturbation modes 
oc exp [— 2Ht] and oc exp [3 (7 — 1)]. While the first mode 
exhibits the same behaviour (exponential damping) as 
the dominating mode for A = and 7 = 0, the behaviour 
of the second mode strongly depends on 7. It decays ex- 
ponentially for any 7 < 1, e.g. for K-matter with 
7 = 2/3, but grows for 7 > 1, e.g. for radiation with 
7 = 4/3. Comoving perturbations in relativistic matter 
in a A-dominated universe are exponentially unstable. A 
behaviour such as this was first found in |2^] by directly 
solving a second-order equation of the type (|l^). A A- 
dominated universe differs from the previous ones with 
A = in so far as the background dynamics according to 
( |l0| ) is completely determined by A in the case considered 
here ( up <C A). On the other hand, differentiating Eq. 
(||) in order to obtain Eq. ([l4]) the cosmological constant 
drops out and does not enter the perturbation dynamics 
explicitly. There are perturbations only in the compo- 
nent which is dynamically negligible in the background 
but not, by definition, in A which, however, determines 
the evolution of a and O in formula ( p3| ) . The mentioned 
instability refers to this kind of perturbations. 
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V. SUMMARY 

We have presented a compact, covariant version of lin- 
ear cosmological perturbation theory for media charac- 
terized by an energy-momentum tensor of the perfect 
fluid type. The general solution of the large-scale pertur- 
bation dynamics in a spatially flat universe was obtained 
as a closed time integral for comoving energy density per- 
turbations, resulting in a comprehensive picture of the 
cosmic mode structure for different cosmological epochs 
including exponential and power-law inflation as well as 
a K-matter period, the standard FLRW cases and matter 
perturbations in a A-dominated universe. 
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